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Abstract. We present observations and analyses of seven 
Galactic O stars of type 06 and earlier. The analyses are 
carried out using NLTE plane-parallel, hydrostatic models as 
well as NLTE spherical models with mass-loss. With detailed 
calculations for the former and simulations for the latter, it is 
shown that the flux blocking due to UV metal lines is important 
for these objects, in agreement with previous studies, and the 
way the mechanism operates is explained. We find that the 
plane-parallel, hydrostatic unblanketed model atmospheres 
have increasing difficulties in fitting the early-type spectra of 
massive stars, and for 50000 K and above a fit seems to be 
impossible. The gravities derived are relatively low even for 
the luminosity class V stars. These objects also show the mass 
discrepancy found in earlier studies, indicating that sphericity 
and mass-loss are important, even at their higher gravities. 
We then perform an analysis using spherical models with 
mass-loss. It is found that gravities should be increased by 
0. 1-0.25 dex, reducing, but not solving, the mass discrepancy. 
We show that spectroscopic masses are in better agreement 
with the theory of radiatively driven winds than evolutionary 
masses are. A helium abundance larger than solar is also 
obtained for most objects. 

Some additional effects (partly related to present appro- 
ximations) that have an influence in our analyses are studied. 
It is found that He II A4200 is less sensitive to details of the 
model calculations than HellA4541 and thus it is preferred 
for temperature determinations, with the consequence of lower 
effective temperatures. It is shown that the fits to He II A4686 
are improved when the upward rates of the Hell resonance 
lines are reduced (with respect to the conventional treatment 
adequate for lines formed in expanding atmospheres), either 
by setting them in detailed balance or by artificially adding 
extra opacity sources that simulate line blocking. The Hell 
blend with H Q is also affected. 
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Some stars of our sample have such high mass-loss rates that 
the derivation of gravities from the wings of Balmer lines, in 
particular H 7 , becomes doubtful. For the most extreme objects, 
the mass-loss rates needed to fit H Q are different from those 
needed to fit H 7 , by a maximum factor of two. 

From the point of view of individual stars, we have ana- 
lysed some of the most massive and luminous stars in the 
Milky Way. According to our analysis, three of them (Cyg OB2 
#7, HD 15 570 and HD 15 558) have particularly large initial 
masses, close to or in excess of 100 M Q . Finally, the least lu- 
minous object in our sample, HD 5 689, could have been erro- 
neously assigned to Cas OB7 and might be a runaway star. 

Key words: Stars: atmospheres - Stars: early-type - Stars: evo- 
lution - Stars: fundamental parameters - Stars: mass-loss 



1. Introduction 

Massive hot stars are relevant objects in many areas of astro- 
physics. They are the main sources of photoionization of their 
surrounding interstellar medium, and they contribute to the 
chemical and dynamical evolution of the host galaxy through 
their stellar winds and supernova explosions. They are also the 
precursors of Wolf-Rayet stars and LBVs and constitute excel- 
lent probes of evolutionary models. For these reasons, advances 
in our understanding of these objects have an impact on other 
fields of astrophysics. 

In spite of all this, few attempts have been made to analyse 
massive O stars quantitatively. The most comprehensive study 
up to now has been that of Herrero et al. (1992 hereafter Paper 
I), who analysed 24 stars from spectral types B0 to 05. De- 
tailed quantitative analyses of massive stars of early spectral 



type have been even more scarce. Conti & Frost (1977) first 



made a systematic analysis of the earliest spectral types, com- 



paring with theoretical predictions by Auer & Mihalas (1972), 



but only for gravities of log g= 4.0. Later, Kudritzki (1980) stu 



died HD 93250, an Q3 V star, using models simila r to those of 
Auer & Mihalas jl972| ), and Kudritzki et al. ( |l983| ) studied the 
spectrum of £ Pup and showed that the star has a low gravity 
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(of the order of log 5= 3.5) and a high helium abundance (e= 
iV(He)/(JV(H)+iV(He))= 0.14). Kudritzki & Hummer dl990|) 



list 15 stars earlier than spectral type 06 (only three classified 
as supergiants), with parameters determined using the same 



methods. Puis et al. (1996) list 22 stars earlier than 06 in the 



Milky Way and the Magellanic Clouds, still taking advantage 
of the optical analysis, but already incorporating the effects of 



sphericity and mass-loss. Pauldrach et al. (1993) analysed Mel 



nic k 42 an d £ Pup, using only the UV spectrum, and Taresch et 
al. ( 1997 ) made a very detailed analysis of HP 93 12 9 A using 
the UV and FUV spectrum. Finally, de Koter ( 1998 ) has used 
the UV Ov line at 1371 A to determine the temperature of 
very hot stars in R136a. 

All these studies have revealed a number of problems in the 
analyses of these stars, related to our incomplete understanding 
of these objects. In Paper I, the so-called helium and mass dis- 
crepancies, already present in the previous literature (see Ku- 



dritzki et al. 1983; Voels et al. 1989; Groenewegen et al. 1989; 



Herrero et al. 1990) were shown to be systematic. These refer 



to the discrepancy in the values of the stellar mass and the pho- 
tospheric helium abundance obtained from the analysis of the 
spectrum using state-of-the-art model atmospheres and evolu- 
tionary models. The explanation of these discrepancies is still 



unclear (for recent working directions, see Howarth 1998) 



In Paper I we already noted the correlation between the 
mass discrepancy and the distance of the star to the Edding- 
ton limit, indicating that the plane-parallel geometry and the 
hydrostatic equilibrium assumption could be the reason for the 
low stellar masses derived. However, the use of wind tech- 
niques in the same work already indicated that the discre- 
pancy could be reduced, but not solved, by including mass-loss 
and sphericity effects. This was later confirmed in an analy- 



sis of HDE 226 868 (Herrero et al. 1995), where the authors 



used Unified Model atmospheres to determine the mass of this 
star, the optical counterpart of Cygnus X— 1, by combining the 
spectroscopic analysis with the orbital data. A similar result 



has been obtained by Israelian et al. (1999) in an analysis of 



HD 188 209. The inclusion of mass-loss and sphericity did not 
seem to hav e any effect on the helium discrepancy either (He- 
rrero et al., 1995 ) in spite of the variations that strong winds 
could introduce in the helium profiles as compar ed to static, 
plane-parallel atmospheres (Schaerer & Schmutz 1994 ). One 
of the reasons for this weak influence was the fact that the study 
in Paper I was limited to spectral types of 05 or later, because 
it was found that above 40 000 K the neutral helium singl et and 
triplet lines gave different stellar parameters. Herrero ( 1994 ) 
suggested that this was due to the neglect of the so-called line- 
blocking, the UV background opacity due to metal lines, during 
the line formation calculations. Also the inclusion of microtur- 
bulence in these calculations can reduce this difference (McEr- 



lean 1998; Smith & Howarth 1998 Villamariz & Herrero, in 



prep.). 

In this paper we study a few stars of early spectral type in 
an attempt to cover several objectives. First, we would like to 
extend our sample from Paper I towards earlier spectral types 
and thus cover the whole region of interest in the HR diagram 



with plane-parallel analyses to see their complete behaviour. 
Although the plane-parallel models will have difficulties in ex- 
plaining even the optical spectrum of these very hot stars (see 
Sect. Q), this first step is needed for the subsequent application 
of more sophisticated models, which will use the experience 
gained and the parameters obtained as input. An analysis of 
line-blocking effects is mandatory here, as it has been seen to 
have an influence for the higher temperatures. 

Then, we will repeat the analyses using spherical models 
with mass-loss, and present the study of some effects that both 
can influence the determination of stellar parameters and that 
will help us to gain new insight into the physics of these stars. 
Having these parameters determined, we can try to establish 
conclusions with respect to the use of different model atmo- 
spheres and techniques for the analysis of very early spectral 
types. We will compare the results of the spherical models with 
mass-loss with those from plane-parallel, hydrostatic model at- 
mospheres as well as with the results obtained using t he som e- 
what approximate technique employed by Puis et al. ( |1996| ). 

In Sect. ^ we present the observations. The effects of line- 
blocking in plane-parallel models are treated in Sect. ^| and the 
spectral description and plane-parallel analysis are considered 
in Sect. fij Sect. § shows the analyses performed with spheri- 
cal models with mass-loss, while Sect. f| contains a qualitative 
study of some effects of interest which explain some difficulties 
found in the preceding section. Then we present our discussion 
(Sect. and conclusions (Sect. fy. 

2. The observations 

The observations were carried out with the 2.5 m Isaac Newton 
Telescope at the Roque de los Muchachos Observatory on La 
Palma during two different observational runs, in 1991 Septem- 
ber and 1992 August. We have observed the spectral region 
between 4000 and 5000 A, and the region around H a . The In- 
termediate Dispersion Spectrograph was used with the H 2400 
B grid in the blue and the H 1800 V grid in the red, attached 
to the 235 mm camera, which resulted in spectral resolutions 
of 0.6 and 0.8 A, respectively, measured on the Cu-Ar arc. The 
main difference between both runs was the size of the CCD 
detector, with a larger wavelength coverage per exposure (400 
A instead of 220 A in the blue) during the second run. Ta- 
ble I gives the stellar identification (usually the HD number), 
star name, OB association to which the star belongs, spectral 
type and the night in which the stellar spectrum was obtained. 

The reduction of the data was done following standard pro- 
cedures. We used both IRAF^and own software developed in 
IDL. The latter was also used for the spectral analysis. The S/N 
ratio of the reduced spectra depends on the spectral range, but 
is usually about 200 in the H 7 region. 



The IRAF package is distributed by the National Optical Astro- 
nomy Observatories, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under contract with the Na- 
tional Science Foundation 
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Table 1. Journal of observations. September dates refer to 1991, and August dates to 1992. The night of 9/10 August was 
dedicated to H Q observations. All spectral classifications are taken from Walborn (1972, 1973), except that of HD 5 689, which 



is taken from Garmany & Vacca ( 1991 1. Properties indicated in the last column have been taken from Mason et al. (1998) 



Order 


Identification 


Spectral Type 


Name 


Assoc. 


Observing nights 


Notes 


1 


Cyg OB2 #7 


03 If* 




Cyg OB2 


6/7, 9/10 Aug 




2 


HD 15 570 


04 If+ 




Cas OB6 


6/7/8, 9/10 Aug 




3 


HD 15 629 


05 V((f)) 




Cas OB6 


28/29 Sep; 7/8, 9/10 Aug 




4 


HD 15 558 


05 111(f) 




Cas OB6 


28/29 Sep; 9/10 Aug 


binary 


5 


HD 14947 


05 If+ 




Per OBI 


28/29 Sep; 9/10 Aug 




6 


HD 210 839 


06 I(n)fp 


ACep 


Cep OB2 


28/29 Sep ; 9/10 Aug 


runaway 


7 


HD 5 689 


06 V 




Cas OB7 


8/9, 9/10 Aug 





3. Line blocking in plane-parallel models 

In Paper I we found that the He I singlet and triplet lines of 
O stars hotter than 40 000 K indicated different stellar param- 
eters when used for the determination of stellar temperatures. 
Following the results of Pauldrach et al. (1993), who stressed 
the importance of line blocking for the wind ionization struc- 
ture and emergent flux, Herrero ( 1 994 > showed that this diffe- 



rence is considerably reduced when UV metal line opacity is 
included. The author attributed this to the different effect of the 
modified UV radiation field on the He I occupation numbers. 
However, no further proof was given there. 

Here we want to discuss in a more detailed way the effect 
of this line blocking, as it is important for the analysis of the 
stars listed in Table [| 

We have included the line list of Pauldrach et al. ( 1993| ) be- 
tween 228 and 912 A. In this region the ionization of He I and 
H takes place. Thus, we can expect an effect on the occupation 
numbers of these two ionization stages. The line list comprises 
roughly 14 000 lines from 26 different elements and 1 37 ioniza- 
tion stages taken into account in the stellar wind NLTE calcula- 
tions by Pauldrach et al. (1993). We calculate their occupation 
numbers in LTE. Although this is a rough approximation, we 
expect that its inclusion will already give us the major part of 
the effect. 

We will use a model at T cS = 40 000 K, log g = 3.40 and e 
= 0.09 to illustrate the effects of line blocking. In Fig. 1 we plot 
the emergent flux between 228 and 912 A. As we can see, a 
considerable fraction of the flux has been blocked. This energy 
will appear at other wavelengths. However, as most of the flux 
escapes between 912 and 2 000 A, this has no real influence 
redward of 912 A. It has to be pointed out that Fig. 1 is only 
illustrative. It is not the emergent flux which is important here, 
but the mean intensity of the radiation field at the depths where 
the He I continuum becomes optically thin. The effect, how- 
ever, is similar. 

The most important effect is produced in the occupation 
numbers of the ground level of He I. As a consequence of 
the flux blocking, the ionization from this level (its ionization 
edge lies at 504 A) is considerably reduced, and its population 
largely increases (see Fig. 2). This does not have any noticeable 
effect on other ionization stages, since He I is very scarcely 



populated. The lower level of the He I AA4387, 4922 singlet 
lines, the 2p 1 P° level, is directly connected to the ground level 
(which also belongs to the singlet system) through a radiative 
transition at 584 A, and thus partially follows their changes and 
also increases its population. On the contrary, the lower level of 
the He I A4471 triplet line (the 2p 3 P° level) is only weakly or 
indirectly coupled to the He I ground level. In addition, the 
ground level of the triplet system (to which 2p 3 P° is strongly 
coupled) is dominated by its ionization and recombination at 
around 2600 A, and it is not affected by the line-blocking. As a 
result, the behaviour of 2p 3 P° follows that of the He I ground 
level to a much lesser extent than does the 2p 1 P° level (as 
do all triplet levels compared to their singlet counterparts). We 
can see this behaviour in Fig. 2, where we have additionally 
marked the formation depths of the center of the He I AA 4387, 
4922 and 4471 lines and the continuum. We can see that the 
populations of these levels increase over the whole region of 
formation, the changes being larger for the 2p 1 P° level. In 
addition, the formation region of the singlet lines is more ex- 
tended than in the case without line-blocking (an effect that is 
not seen in the figure). The combination of larger changes over 
a larger line formation region produces the stronger variations 
of the singlet lines compared to the triplet one. 

Figs. 3 and 4 show the variations in the line profiles of 
He I AA4387, 4471. We see that the changes are much less 
in the second one, and thus this line should be preferred for 
analyses at high temperatures. At lower temperatures and grav- 
ities, however, and due to the so-called dilution effect in this 
line, the use of the singlet lines is preferable, even if the calcu- 
lations are made without line-blocking. This effect refers to the 
fact that the fit of the He I A4471 line becomes worse when 
going from dwarfs to supergiants, while the rest of the lines re- 
tain a good fit quality (see Voels et al. 1989). Although still not 
completely clarified, Smith & Howarth (1998 ) recently claimed 
that microturbulence could be the cause. 

At lower temperatures or larger gravities, collisions play a 
stronger role, and the effects of line blocking are lower. This ex- 
plains the behaviour of the corrections to the stellar parameters 



found by Herrero (1994 especially Fig. 3) when including line- 



blocking, and basically consist of obtaining lower temperatures 
when line-blocking is not taken into account. The amount of the 
correction will depend on the model temperature at a given gra- 
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Fig. 1. The Eddington flux in the region between 228 and 1200 A, with and without approximate line blocking, for the model at 
T eff = 40 000 K, log g = 3.40 and e = 0.09. 
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Fig. 2. The changes in the occupation numbers of the ground level of He I (full line), the lower level of the singlet lines 
He I AA4387, 4922 (2p dashed line) and the ground level of the triplet line He I A4471 (2p 3 P°, dashed-dotted), and 
the formation depths of the line centers and the continuum. Negative numbers mean that the atomic level populations calculated 
with line-blocking are larger. The peak near log m = -2.5 is a numerical artefact in the convergence of the model without line 
blocking and does not appear in other models. The model parameters are the same as in Fig. 1. 

vity, being larger for higher temperatures. It has to be stressed, our typical error boxes (±0.1 in log g and ±0.03 in helium 
however, that the corrections are large only if we use the singlet abundance) have been adopted in the course of the parameter 
lines for the temperature determination without line-blocking, determinations. 
Using the triplet line results in a smaller temperature than us- 
ing the singlet ones (by 500-1 000 K at temperatures around 
45 000 K). Gravity and helium abundance do not significantly 
vary, although sometimes variations of the parameters within 
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Fig. 3. The helium line He I A4387 with and without line- 
blocking (solid and dashed lines, respectively) for the same 
model parameters as in Fig. 1. He I A4922 behaves the same 
way. 
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Fig. 4. As Fig. 3, but for the He I A4471 line. 



4. Spectral analysis using plane-parallel, hydrostatic 
models 



Fig. 5 shows the spectra 
5000 A, whereas Fig. 6 
Before performing 
model atmospheres, we 
velocity displacements, 
crucial in the very early 
commonly used for the 



of all observed stars between 4000 and 
shows the spectra around H a . 
the comparison with the theoretical 
have corrected the spectra from radial 
Correction for the radial velocities is 
type stars, because the H and He lines 
spectral analysis may be contaminated 



by wind effects that fill their red wings. However, these cor- 
rections are particularly difficult: the cores of strong lines may 
be affected by the wind and only a few weak metal lines are 
present. Rotation adds a new handicap, as it broadens lines, 
making them shallower and favouring blends. In addition the 
limited spectral range of each single frame limits the number 
of suitable lines on individual spectra. Thus the first difficulty 
is to find a set of lines appropriate for the measurement of the 
radial velocity. 

We have discarded all H and He lines. This is already ne- 
cessary, since sometimes we can see a trend in the measured 
radial velocity with the excitation potential of the line, indi- 
cating that the line cores are formed in higher layers moving 
faster as the excitation potential decreases. 

Thus the radial velocity correction is based on lines of Si 
(Si iv and Sim), N (Nv, Niv and Nm), O (On) and C 
(C III). Sometimes, a few lines are in emission and not in ab- 
sorption. These lines are used only if they give results concor- 
dant with the rest of the metal lines (usually, this is the case). 
The typical uncertainties we obtain are ±20 km s _1 , which is 
about half of the spectral resolution. This large error is due to 
the existence of only a few lines, in addition with broad cores, 
and reflects the dispersion in the individual values. Within these 
limits, we have sometimes displaced the fitted lines (also in the 
analysis with spherical models), when it was clear that this was 
the only way to obtain a good fit. Table 2 lists the radial velo- 
cities in the kinematical Local Standard of Rest. 

We then determined rotational velocities following the 
same procedure as indicat ed in Paper I. Our values comp are 
well with those by Penny ( |1996[ ) and Howarth et al. ( |1997| ) for 
objects in common (see discussion about individual objects). 
Rotational velocities are given in Table 2 together with the pa- 
rameters determined for each star. In this table, temperatures 
are given in thousands of Kelvin; gravities are corrected for 
centrifugal force effects (see individual discussions for model 
parameters); e is the helium abundance by number with res- 
pect to the total number of H and He atoms; V is the inte- 
gral of the stellar flux over A, weighted by the V-filter function 



of Matthews & Sandage (1963), used to calculate stellar radii 
from the model atmospheres (see Kudritzki 1980| , or Paper I); 
My is the absolute magnitude. This has been obtained using 
the photometry and colours from Hipparcos for HD 15 570, 
HD 15 629,HD 15 558, HD 14947 and HD 210 839 combined 
with th e extinction laws and distances from Garmany & Stencel 
( |1992| ). For Cyg O B2 #7, M v has been taken from Massey & 
Thompson ( 1991 1, and for HD 5 689 has been taken from Gar- 
many & Stencel ( 1992 ). The evolutionary masses hav e bee n 
obtained from the evolutionary tracks by Schaller et al. ( 1992 ). 

Except for the line-blocking, the models are the same ones 
used in Paper I. These are H/He, , plane-parallel models, in 
hydrostatic and radiative equilibrium. The line fit is also made 
in the same way. We first fit H 7 , obtaining the best gravity at 
a given T e g, The same procedure is followed for the helium 
lines, at a normal He abundance (i.e., 0.09 by number). The lo- 
cus where all lines intersect in the log T c ff- log g diagram is 
taken as giving the stellar parameters. If the lines do not inter- 
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Fig. 5. The blue stellar spectra. The relative fluxes have been arbitrarily displaced in steps of 0.5 for the sake of clarity. Wave- 
lengths are in A. 



sect, the helium abundance is changed. The helium abundance 
giving the smallest intersection region is adopted as the one ap- 
propriate for the star. In previous studies, the helium lines used 
were He II AA 4541, 4200 and He I AA 4387, 4922, whereas 
He I A 4471 was used only for dwarfs or high-temperature 
stars. As this is the case in the present study, this is the line 
we use here. In addition, for the reasons explained before, we 
give it a larger weight than for the singlet lines. 

The errors are similar to those quoted in Paper I, with 
±1500 K in T cB , ±0.1 in logg and ±0.03 in e. This pro- 
duces errors of ±0.06, ±0.19 and ±0.22 in \og(R/R Q ), 
\og(M / Mq) and log(L/X Q ), respectively, when adopting an 
uncertainty of ±0.3 mag for My. 

We now describe the line fits of each star independently, 
and the spectral features, if appropriate. In Fig. 7 we show 
the theoretical HR diagram for these stars, with the parameters 
listed in Table 2. 
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4.1. HD5689, 06 V 

This is the only star for which we can clearly see lines of He I 
other than He I A4471. Its line fit, shown in Fig. 8, is obtained 
at T eff = 40000 K, log 3 = 3.40, e= 0.25 and V T sini= 250 km 



Fig. 7. The programme stars on the Hertzsprung-Russell di- 
agram after the plane-parallel analysis. The meanings of the 
symbols are those indicated in the lower left corner, together 
with typical error bars. The theoretical tracks are from Schaller 



The rotational velocity, however, is uncertain, as in this et al. (1992). The numbers to the left of the zero-age main 



cases we only have the He lines (especially He I). The given 
value can be seen as an upper limit to the projected rotational 
velocity. We must correct for centrifugal force by adding the 
term (V r sin i) 2 /R, since the measured value is the effective 
gravity reduced by the centrifugal acceleration, which has to 



sequence indicate the initial stellar masses in units of solar 
masses. L is in solar units. 
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Fig. 6. As Fig. 5, however for the red spectra 

Table 2. Parameters determined for the programme stars using plane-parallel models. Note that the parameters of Cyg OB2 #7 
and HD 15 570 are only indicative, and could not be determined with the plane parallel models. Temperatures are in thousands of 
Kelvin. M s , M cv and Mo are respectively, the present spectroscopic and evolutionary masses, and the initial evolutionary mass, 
in solar units. The last column indicates whether we have formally a mass discrepancy considering an error of 0.22 in log(M s ). 
V x sini and w ra( j are given in km s _1 . Semicolons in v xa d indicate that the uncertainty in this value is larger than ±20 kms -1 



Star 


"rad 


T c ff 




e 


V r sini 


V 


M v 


R/Rq 


log(L/L ) 




M ev 


Mo 


M discr.? 


Cyg OB2 #7 


+35 


51.0 


3.66 


0.12 


105 


-29.655 


-6.20 


16.7 


6.23 


46.2 


111.8 


114 


Yes 


HD 15 570 


-20 


50.0 


3.51 


0.15 


105 


-29.553 


-6.73 


22.0 


6.44 


51.1 


139.1 


142 


Yes 


HD 15 629 


-50 


48.0 


3.81 


0.09 


90 


-29.571 


-5.52 


12.7 


5.89 


37.5 


69.9 


71 


Yes 


HD 15 558 


-35 


46.5 


3.71 


0.07 


120 


-29.512 


-6.28 


18.5 


6.16 


63.8 


91.7 


97 


No 


HD 14947 


-35: 


45.0 


3.53 


0.15 


140 


-29.415 


-5.69 


14.8 


5.91 


26.5 


65.8 


68 


Yes 


HD 210 839 


-70 


41.5 


3.47 


0.25 


250 


-29.357 


-6.17 


18.9 


5.98 


38.7 


67.0 


71 


Yes 


HD 5 689 


-65 


40.0 


3.57 


0.25 


250 


-29.337 


-4.19 


7.7 


5.13 


7.8 


30.1 


31 


Yes 



be added here in order to derive the masses. Corrected for cen- 
trifugal forces the gravity in Table 2 is log g = 3.57. As usual for 
fast rotators, the helium abundance obtained is very large. The 
radius, however, is very small, and the derived spectroscopic 
mass is much lower than the predicted evolutionary one. These 
facts reflect certain problems in the set of parameters for this 
star. T he star has been classified as 06 V by Garmany & Vacca 
(1991). This classification is not in complete a greem ent with 
the star belongin g to C as OB7 (see Humphreys 1978 , or Gar- 
many & Stencel 1992), since the absolute magnitude derived 
from the association distance corresponds to a less luminous 
object than an 06 V star by one magnitude, which causes the 
small radius derived (see Table 5 of Vacca et al. 1996| ). The si- 
tuation is of course much worse if, based on the low gravity, 



we associate this star with a giant and adopt parameters cha- 
racteristics of a luminosity class III object. Thus we have two 
possibilities: 1) the star is a main sequence star that belongs to 
Cas OB7, with a smaller radius than usual for its spectral clas- 
sification, in spite of the large rotational velocity; in this case, 
the low gravity is difficult to explain, or the models are giving 
us completely wrong parameters for this star; or 2) the star has 
parameters typical for a luminosity class III object. In this case, 
the star would probably not belong to Cas OB7, although we 
have found no indications in the literature about this possibility. 
Other combinations of the above arguments are also possible, 
but the question that something is non-standard with the star 
remains. This case is very similar to the one we had in Paper I 
for HD 24912, which is a known runaway star. Using the Oort 
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4540 4545 
igth (Angstrom) 



Fig. 8. The fit to the spectral lines in HD 5 689 using the plane- 
parallel, hydrostatic models, with the parameters given in Table 
2. From top to bottom and left to right: H 7 , Up, He I AA4387, 
4922, He II AA4200, 4541. 



constants given by Lang (1980 1, we obtain a peculiar velocity 
of 39 km s _1 for HD 5689, which, given our uncertainties, does 
not allow us to decide clearly whether it is a runaway or not, 
if we adopt the conventional limit of 30 km s _1 for runaway 
stars. 



4.2. HD 210839, 06 1(n)fp 

This is A Cep, a well known fast rotator. The line fit (see Fig. 
9) is obtained at T cS = 41500, log 5 = 3.40, e = 0.25 and 
V T sini= 250 km s _1 . The rotational velocity is again uncer- 
tain and could be lower, as we suggest later. This is also indi- 
cated by the values found by Penny (1996; 214 km s _1 ) and 
219 km s _1 ). With the centrifugal force 



Howarthet al. (1997 



correction, logy increases to 3.47 (see Table 2). The line fit is 
the best agreement we could find between all the He lines. The 
parameters are very similar to those of HD 5 689, (except that 
now the mass discrepancy is only the usual factor 2) but the 
spectrum shows important differences. A Cep displays Of fea- 
tures, and has strong emission in H Q . This has to be attributed 
to the large difference in luminosity (see Fig. 7 and Table 2). In 




4585 438? 4584 4385 4588 4550 459? 4394 4460 4465 4470 4475 4480 4485 




0.80 .... i .... i .... i .... i .... i .... I n 75 
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4535 4540 4545 
Wavelength (Angstrom) 



Fig. 9. As Fig. 8, for HD 210 839. Here we show He I A 447 1 
instead of He I A 4922. 



Fig. 7, HD 210839 occupies the position of an evolved lumi- 
nous star, already within the instability strip predicted by Kiri- 



akidis et al. (1993 1. In agreement with this location, A Cep is 



known to show profile variations due to non-radial pulsations, 
which however are below the accuracy of our optical spectra 
(cf. Fullerton et al. |1996|, de Jong et al. |1999[). (Note, that we 



have seen temporal variations in the U a profile which was ob- 
served several times. In contrast, however, two spectra taken in 
the region from 4200 to 4600 A showed no significant varia- 
tions in the lines used for the spectral analysis.) Finally, let us 
remark that the IR spec trum o f A Cep has been analysed re- 
cently by Najarro et al. (1998 1, who obtained parameters very 
similar to ours (except in the temperature, for which they find 
a value lower by 4000 K, see discussion). 

4.3. HD 14947, (95//+ 

The spectral lines of this star are fitted at T c s = 45 000 K, log g 
= 3.50, e = 0.15 and V r sini = 140 km s _1 , which agrees with 



the value of 133 km s 1 given by Penny (1996) and Howarth 



et al. (1997). The spectral line fits are shown in Fig. 10. The 



star is an e xtreme Of (without being a transition object, see 
Conti et al. |l995| ), with a large emission in Nlll A4630-40 
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and Hell A4686, and also in H Q . We can also begin to see 
Si IV A41 16 in emission and the N V AA4604, 4620 lines in ab- 
sorption. The line fit, shown in Fig. 10, can be considered as 
acceptable. The fit of the He I A4922 line indicates that the 
predicted line is a little too strong as compared to the observa- 
tions, but the difference is very small compared to variations of 
the line within the error box. The He abundance is not as large 
as for the two previous objects, but is still larger than solar. The 
spectroscopic mass is more than a factor of two smaller than 
the evolutionary one. 

4.4. HD15558, 05111(f) 

This star is a binary, but we expect the spectrum not to be conta- 
minated, as it is a single component in a well separated system 



(Mason et al. 1998). The best line fit, shown in Fig. 11, is ob- 
tained at T off = 46500 K, log 5 = 3.70, e = 0.07 and V T sini = 
120 km s -1 . The rotational velocity value agrees with the 123 



km s of Howarth et al. (1997) although it departs slightly 



from the one given by Penny (1996) of 147 km s . The Of 



features are weaker than in the previous star, as is the emission 
in Si IV A41 16, but note that the gravity is now larger. The line 
fit with hydrostatic models is very difficult. Only the wings of 



H 7 and Hp are fitted, together with He II A4541. The singlet 
He I lines are again too weak and noisy, but in this case also 
He I A447 1 cannot be fitted. We cannot simply attribute the dis- 
tortion in the blue wing of He I A447 1 to binarity. If this were 
the case, we should see the singlet lines too (if the companion 
is relatively cool) or distortions in the He II lines (if the com- 
panion is relatively hot). The lack of a good fit in any He I line 
makes the parameter determination much less accurate. From 
the plane-parallel hydrostatic models this is the most massive 
star, and in fact the mass discrepancy is comparatively low, the 
difference between the spectroscopic and evolutionary masses 
being only 30%. 

4.5. HD 15629, 05 V((f)) 

The best line fit for this star is obtained at T cS = 48 000 K, 
log g = 3.80, e = 0.09, with V x sin i = 90 km s _1 . This last value 
agrees within ±1 km s _1 with the values quoted by Penny 



(1996) and Howarth et al. (1997). In spite of the high tempera 



ture, the larger gravity allows us the calculation of hydrostatic 
models. The line fit is shown in Fig. 12. We can see that the 
fit is good for the wings of H 7 and H^, and for He I A4471 
and Hell A4541, but it is bad for He I A4387, and also for 
He I A4922, that is not shown in Fig. 12. We have given much 
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less weight to these two lines, however, because they are very 
weak and noisy. In the case of He I A4922, the local conti- 
nuum has also been placed too low. Correction of this would 
bring the calculated line into agreement with the observations. 
More worrying is the lack of fit in the He II A4200 line. The fit 
of this line is bad, and comparable to that of the He II A4686 
line, for which we expect a bad fit with hydrostatic models. 
For this reason we preferred to give more weight to the fit of 
He II A4541 than to He II A4200. In spite of the large gravity 
and the luminosity class V, this star also shows the mass discre- 
pancy, which was not the case for less luminous stars in Paper 
I, where we found no significant mass discrepancy for stars of 
high gravities. This indicates the increasing role of radiation 
pressure. 

It is interesting to note the case of HD 15 629, HD 14947 
and HD 210 839. Within our error bars, these stars could rep- 
resent different evolutionary stages of a star of initially around 
70 M Q following standard evolutionary tracks with mass loss 
but without rotation. Within this scenario, it is impossible to 
explain the higher He abundances of the two cooler stars. This 
can be an indication that rotation plays a strong role in stellar 
evolution, since there is no other known mechanism that might 
substantially modify the atmospheric He content of a single star 
like HD 15 629 in only two Myr. (Another possible scenario for 
an increased He abundance, close binary interaction, can most 
probably be discarded in the case of these three particular ob- 



jects, see Mason et al., 1998). Turning the argument around, 



if rotation does not play a significant role in the evolution of 
single massive stars, HD 15 629 could become rather similar to 
first HD 14947 and then HD 210839 in one to two Myr. Un- 
fortunately, we cannot wait to confirm this hypothesis. 

4.6. HD15570, 04 If+ 

This star is hotter than the former one, and its Of signatures are 
str onger (again, without being a transition object, see Conti et 
al. 1995), indicating a lower gravity. This is confirmed by ins- 



pection of the Balmer series, especially the H a profile, which 
is strong in emission. From 4600 to 4750 A we see a broad 
emission feature with emission lines typical of Of stars, and 
both Si IV lines neighbouring Hg are in emission, indicating a 
large luminosity, which is confirmed. There is also absorption 
in N V AA4604, 4620. It is by far the most luminous star in the 
sample, and we were not able to fit the spectrum properly for 
this object. The extreme character of the features allows only a 
crude guess of the stellar parameters, and we do not show any 
line fit. Apart from the projected rotational velocity of 105 km 
s , we can only say that the temperature is close to 50000 K 
and the gravity should be of the order of log g ~ 3.50 or even 
less. With the given luminosity, the evolutionary mass is very 
large, indicating an initial mass in excess of 140 M Q , which 
would make HD 15 570 one of the most massive and luminous 
stars known in the Milky Way. The spectroscopic mass is much 
lower, although we have to consider that the gravity has only 
been "guessed" so far (see next section). There is a chance that 
the star is helium enhanced, and we adopt e = 0.15. 
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Fig. 12. As Fig. 9, for HD 15 629. 
4.7. CygOB2#7, 03 If 

This is a very hot star. He I A4471 is only marginally seen, and 
the emission in N IV A4058 and the absorptions in N V AA4604, 
4620 are strong. We again find the broad emission feature with 
emission lines typical for Of stars, and both Si IV lines neigh- 
bouring Hs in emission. In spite of our efforts, we were unable 
to fit the star with plane-parallel models, the required gravi- 
ties being too low. As for HD 15 570, we had to extrapolate the 
plane-parallel parameters, and we do not show a line fit. The 
final parameters we adopt are T c s = 51 000 K, \ogg = 3.65, e 
= 0.12 and V r sini = 105 km s -1 . We did not succeed in calcu- 
lating models in this range (log g had to be larger at least by 
0.1), and spherical models including mass-loss are obviously 
required. This star is also very luminous and massive, with an 
initial mass in excess of 1 10 Mq, but with a large mass discre- 
pancy. 

5. Analysis with spherical models with mass-loss 

The difficulties in the analysis of some of the spectra discussed 
in the previous section result from the neglect of sphericity 
and mass-loss effects. Thus, we d ecided to use the program 
described by Santolaya-Rey et al. ( 1997 > to account for these 
effects. Briefly, this program constructs a unified model of the 
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stellar photosphere and wind regions, with a number of stan- 
dard (stationarity, homogeneity) and non-standard assumptions 
and approximations (the most important of the latter concerns 
the way the temperature structure is obtained, which results in 
well approximated, but not forced, flux conservation through- 
out the atmosphere and a constant temperature in the wind). It 
solves the line-formation problem in an expanding atmosphere 
with spherical geometry, and Stark broadening is included in 
the final formal solution. 

The atomic models used are the H and He models used 
also in our plane-parallel analysis, with minor changes to adapt 
them to the new program. Line-blocking has not been included 
so far. Its influence, however, will be considered in the next sec- 
tion (at least on a qualitative basis), and has affected the results 
decribed here. 

We have adopted the following procedure. We begin with 
the parameters given by the plane-parallel models (see Table 
2). The values of the terminal velocities are taken from Puis et 
al. (1996), except for HD 15 570 (from Lamers & Leitherer 
19930 and for Cyg OB 2 #7 and HD 5689, where we have 
adopted the Voo - spectral classification relation provided by 
Haser (1995). The last object has been assigned a Voo corres- 



6.5 ■ 



ponding to an 06 III star, which is a compromise between the 



luminosity class V given by Garmany & Stencel (1992) and 



the low gravity obtained from our analysis (see the discussion 
in the preceding section). Note, that this compromise is reason- 
able, since the difference between the mean Voo of 06V and 
061 stars is only 450 km s _1 . 

Also needed are the values of the exponent in the f3— 
velocity law. For HD 15 558, HD 15 629 and HD 210839 we 



adopted the values given by Puis et al. (|1996J). For HD 14947, 
also analysed by Puis et al., we adopted j3 = 1.15 instead of 
1.0 as quoted by Puis et al., for reasons explained later in the 
discussion of HD 14947. For Cyg OB2 #7, HD 15 570 and 
HD 5 689 which were not analysed by Puis et al we used j3= 
0.8, 1.0 and 1.0, respectively. 

For each star we have constructed a small model grid by 
varying gravity, temperature and mass-loss rate. We then tried 
tofitH 7 ,H a and the ratio He II A4200toHel A4471, which 
is our new preferred temperature indicator, following the re- 
sults described in the next section, where we will show that this 
He II line is less sensitive to model details than He II A4541 . Of 
course, the use of a new temperature indicator introduces also 
an additional difference to the results from our plane-parallel 
considerations. 

All other parameters remain fixed at the beginning at their 
values determined from the plane-parallel analysis. When ne- 
cessary, the helium abundance has been changed later. The ra- 
dius merits additional comments. Its value is derived from the 
emergent flux, which has been obtained for a given set of pa- 
rameters (those from the plane-parallel analysis). Changing the 
parameters will change the derived radii, and thus an iteration 
process might be required. We have recalculated the radii of 
our objects with the new parameters, and have obtained an av- 
erage change of only 2% (the maximum change being of only 
4.5%, from 22 to 23 Rq). Thus we used the plane-parallel va- 




Fig. 13. The programme stars on the Hertzsprung-Russell dia- 
gram after our analysis with spherical models with mass-loss. 
Compare the new diagram with Fig. 7. Details are as in Fig. 7. 



lues without further iteration. The finally adopted and derived 
parameters are summarized in Ta ble 3. For those stars which 
have been analysed by Puis et al. ( 1996 ), we give their results 
in Table 4 for comparison. The new position of the stars in the 
HR diagram is shown in Fig. 13. 



5.1. HD5689, 06 V 

We can see the line fit to this star in Fig. 14, for a model with 
T e g = 37 000 K, logg = 3.45 (which increases up to 3.57 due 
to the centrifugal correction) and log M = -6.80. The tempe- 
rature is lower than in the plane-parallel case, in part due to 
the change of our temperature criterion. The gravity, however, 
remains the same (which actually implies an increase with res- 
pect to plane-parallel models, since a lower temperature usu- 
ally demands a lower gravity to fit a given line). Thus, the large 
mass discrepancy found above (see Table 2) cannot be reduced 
by the spherical models, since it originates from the low value 
of the radius, rather than from the low gravity. Neither can we 
decide whether the star has a larger radius in parallel with a 
higher mass-loss rate, since the H Q wind-emission depends 
only on the ratio (M /i? 3 / 2 ), whereas its photospheric compo- 
nent is independent of both. Assuming an absolute magnitude 
of a typical 06 III star (My = -5.78 mag, see Table 6 of Vacca 
et al. 1996| ) results in a radius of 15.9 Rq, but does not change 
the fit quality (it affects only slightly the fit of He I A4471). 
Of course, mass, luminosity and mass-loss rate would be af- 
fected and we would obtain a spectroscopic mass of 30.0 Mq, 
an evolutionary mass of 41 .9 M Q (thus reducing the mass dis- 
crepancy to 30%, a typical value), a luminosity of 5.6, closer 
to the other objects (see also the discussion in Sect, and a 
logarithmic mass-loss rate of -6.17. However, this possibility 
would also imply that the wind in HD 5689 is much less effi- 
cient than in HD 210 839, as both stars would have roughly si- 
milar parameters, except for M. Unfortunately, without know- 
ing the distance to HD 5689 more accurately, we cannot derive 
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Table 3. Parameters determined for the programme stars starting with the parameters from Table 2 and using the spherical models 
with mass-loss. Temperatures are in thousands of Kelvin, velocities in km s _1 and mass-losses in solar masses per year. MWM 
stands for modified wind momentum, and given is log (M VooR ' 5 ), with M in solar masses per year, Vqq in km s _1 and R in 
solar radii. M s , M ev and Mg are, respectively, the present spectroscopic and evolutionary masses, and the initial evolutionary 
mass, in solar units. The last column indicates whether we formally have a mass discrepancy considering an error of 0.22 in 
log(M s ). 



Star 


T c s 


logs 


e 


Rj Rq 


log(L/L Q ) 






log M 


M B 


M ov 


Mo 


MWM 


M discr.? 


Cyg OB2 #7 


50.0 


3.72 


0.18 


16.7 


6.20 


2900 


0.80 


-4.95 


53.0 


104.5 


107 


29.93 


Yes 


HD 15 570 


42.0 


3.80 


0.15 


22.0 


6.14 


2600 


1.05 


-4.75 


112.6 


79.6 


89 


30.14 


No 


HD 15 629 


46.0 


3.81 


0.09 


12.7 


5.81 


3000 


1.00 


-6.13 


37.5 


61.4 


63 


28.70 


No 


HD 15 558 


46.5 


3.86 


0.07 


18.5 


6.16 


2800 


0.75 


-5.40 


89.5 


91.7 


97 


29.48 


No 


HD 14947 


40.0 


3.67 


0.20 


14.8 


5.70 


2400 


1.15 


-5.25 


36.8 


47.5 


49 


29.52 


No 


HD 210 839 


37.0 


3.55 


0.25 


18.9 


5.78 


2250 


0.90 


-5.17 


45.8 


49.3 


52 


29.62 


No 


HD 5 689 


37.0 


3.57 


0.25 


7.7 


5.00 


2500 


1.00 


-6.80 


7.9 


25.4 


26 


27.85 


Yes 



Table 4. Parameters as determined by Puis et al. (1996) for some of the stars in Table 3. Units are as in Table 3. 



Star 


T c tt 


logs 


e 


Rj Rq 


log(L/L ) 


\4o 


P 


log M 


log (M VvoR' *) 


HD 15 629 


47.0 


3.90 


0.07 


14.2 


5.95 


3000 


1.00 


-6.12 


28.73 


HD 15 558 


48.0 


3.85 


0.07 


21.8 


6.36 


2800 


0.75 


-5.14 


29.78 


HD 14947 


43.5 


3.50 


0.15 


16.1 


5.93 


2350 


1.00 


-5.12 


29.65 


HD 210 839 


38.0 


3.65 


0.09 


19.0 


5.83 


2250 


0.90 


-5.28 


29.52 



stronger conclusions. We also see that He II A4541 does not fit 
completely, and in particular the fit of He II A4686 is poor. This 
line always shows a poor fit, with an observed absorption that 
is much stronger than predicted, especially in the blue wing. 
As explained in the following section, this seems to be related 
to lack of line -blocking in our models, and He II A4541 seems 
also to be affected, although to a lesser extent. 

5.2. HD 210839, 06 1(n)fp 

The line fit of this star, shown in Fig. 15, reveals two problems. 
First, the form of H 7 suggest that the adopted rotational veloc- 
ity is too large. A value of 200 km s _1 results in a much better 

The second problem is that we were unable to fit the P 
Cygni form of the H Q -profile. As was pointed out previously, 
A Cep is a possible non-radial pulsator (Fullerton et al. 1996, 



de Jong et al. 1999), which might induce deviations from ho- 
mogeneity (by exciting the line-driven wind instability already 
in the lower wind part, see Feldmeier et al. 1997 ), and the large 
rotation rate might have an additional influence on the wind 
structure and the resulting profile (cf. Owocki et al. 1998 and 



centrate on the red wing of H Q , since the theoretical simula- 
tions result in a blue wing affected by an extra He II absorption, 
which is inadequately described in our present models (see next 
section), and any wind variability becomes much more visible 
in the blue wing, compared to the red one. Again, the predicted 
emission in He II A4686 is much stronger than observed; also, 
the temperature is lower than in the plane-parallel case. How- 
ever, the gravity is larger, and we do not find a mass discre- 
pancy for this star. 



A Cep has been also analysed by Puis et al. (1996), who 



adopted slightly different parameters. The main difference is 
the helium abundance (e = 0.09 instead of 0.25), which has a 
small influence on the gravity, however produces a larger wind 
momentum with a lower luminosity in our present results. (The 
value for e quoted by Puis et al. was not derived from a con- 
sistent spectral analysis, however taken from the literature). It 
is interesting to note that the high abundance favoured by our 
findings is in agreement with Blaauw 's sug gestion that all run- 
away stars are He enriched (Blaauw, 1993 ). Note also that the 
possible runaway nature of HD 5 689 fits within this scenario. 

Finally, we point out that our mass-loss rate coincides well 
with other determinations from H a by Lamers & Leitherer 



references therein; Petrenz & Puis 1996). We decided to con- 1993, but is larger than that derived from radio fluxes (Lamers 



a misprint in Table 8 of Puis et al. (1996) gives a velocity of 100 



& Leitherer 1993) by a factor of 3 



km s 1 for this star, whereas 200 km s was the value adopted in that 
work , and wo uld ag ree with previou s dete rminations (Conti & Ebbets, 



1977, Penny, 



199€ , Howarth et al.. 



1997). A change in the rotational 



velocity would not change the adopted parameters, only the quality 
of the line fit. We show the fit with the value adopted for the plane- 
parallel analysis, although a better fit is possible with the lower V t sin i 
value. 



5.3. HD 14947, 05 If+ 

The analysis of HD 14947 has been hindered by an inconsis- 
tency in the radial velocities derived from the blue and the red 
spectrum (the spectra were taken in September 1991 and Au- 
gust 1992, see Table 1). While in the blue spectrum there are 
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Fig. 14. Fit to the H and He lines in HD 5 689 using a spherical 
model with the parameters given in Table 3. From left to right 
and top to bottom we show H 7 , H Q , He I 4471, Hell 4200, 
He II 4541 and He II 4686. Wavelengths (along the abscissa) 
are given in A. The ordinates give the relative fluxes. 



several metal lines to derive the radial velocity correction, in 
the red spectrum we only have three lines of He II and they 
give a result that is inco mpatible with that of the blue metal 
lines. Mason et al. ( 1998 ) list this star as having a constant ra- 
dial velocity and thus we don't have a clear explanation for the 
radial velocity change. Therefore, the radial velocity correc- 
tion for this star is particularly inaccurate. We have adopted the 
correction indicated by the blue metal lines, accounting for the 
different rest frames as a function of observation date. 



The line fit of H 7 is comparatively poor for this star (see 
Fig. 16). This reflects a problem that we have found for the 
first time, an inconsistency between the fit for H 7 and H a . 
The fit for H 7 would need a mass-loss rate which is only half 
that of H Q . This inconsistency is only weakly dependent on the 
gravity. 

The new temperature derived for HD 14947 is much lower 
than that obtained from plane-parallel models. This is due 
to sphericity together with the new criterion of using He II A 
4200. Again, in agreement with the large mass-loss rate, 



Fig. 15. As Fig. 14, but for HD 210 839. The fit of H 7 improves 
considerably if the adopted rotational velocity is lowered from 
250 to 200 km s" 1 . 



He II AA4541, 4686 are poorly fitted, (see discussion in the next 
section). 

The spectroscopic mass derived is clearly lower than the 
evolutionary one, alth ough t hey agree within the errors. The 
analysis by Puis et al. (1996) gave slightly different stellar pa- 
rameters, the most important difference being the now lower 
temperature implying a lower luminosity by 0.23 dex. 

There is an upper limit of 10~ 4 75 M /yr for the mass- 
loss rate of HD 14947 derived from radio fluxes (Lamers & 



Leitherer 1993). The value we find here is in agreement with 



this upper limit. 

5.4. HD 15558, 05111(f) 

The line fit to HD 15 558 obtained for the parameters given in 
Table 3 can be seen in Fig. 17. The mass-loss rate is rather 
large, in agreement with the high luminosity, especially when 
compared with HD 15 629, a star that could be considered si- 
milar at first inspection since both H Q profiles are in absorp- 
tion. However, the value we obtain here for M is lower than 



that obtained by Puis et al. ( 1996 ) due to the change in the ste- 
llar parameters. The gravity is largely increased with respect to 
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Fig. 16. The line fit to HD 14947 using spherical models with 
mass-loss. We see the comparatively poor fit for H 7 in the core, 
which can be largely improved by reducing the mass-loss rate 
by a factor of two. 



the value derived with plane-parallel, hydrostatic models (by 
0.15 dex), and in fact the derived spectroscopic mass is nearly 
equal to the evolutionary one, thus making HD 15 558 one of 
the few cases for which we do not find any mass discrepancy. 
The generally good agreement in the line fit to HD 15 558 is 
again broken by He II A4686. 

5.5. HD 15629, 05 V((f)) 

The line fit of HD 15 629 is obtained for the parameters given 
in Table 3 and is shown in Fig. 18. We see that the mass- 
loss rate is relativel y low and agrees very well with the value 
given by Puis et al. ( 1996 ). The gravity we obtain is similar to 
that in the plane-parallel case, and thus the new mass is again 
lower than the evolutionary one, although there is formal agree- 
ment considering the error bars. The worst fit cor responds to 
He II A4686. A comparison with Puis et al. ( |1996| ) shows that 
both sets of values are compatible, although the general trends 
(lower temperatures, radii and luminosities in our case) remain. 



Fig. 17. As Fig. 14, but now showing the line fit to HD 15 558. 
5.6. HD 15570, 04 If+ 

This is an extreme object, and the adopted parameters are ac- 
tually a compromise, reflected in very large uncertainties. For 
this star we derive T cff = 42000 K, log 5 = 3.80, e = 0.15, log 
M = -4.75 (see Table 3 for the rest of the parameters). The very 
high mass-loss rate results in an extreme insensitivity of the 
wings of H 7 to gravity variations, as can be seen in Fig. 19, 
where we have plotted the profiles for \ogg = 4.05 and log g = 
3.55. These can be considered as limits beyond which the fits 
begin to become poorer than for the range [3.55,4.05]. Thus we 
have adopted log g = 3.80±0.25. However, it should be borne 
in mind that even such a large error is optimistic, as the uncer- 
tainty of the radial velocity correction has not been included. 
Also the He profiles are largely insensitive to log g variations, 
but they allow us to restrict T c s and e. 

The line fit to this star is shown in Fig. 20. A detailed fit of 
the strong H Q emission of HD 15 570 was not possible, since 
for the given parameters the simulations always showed a dou- 
ble peak due to the extra emission coming from the He II blend 
(see the discussion in the following section) not seen in the ob- 
servations. Thus the adopted criterion was to fit the red wing 
of H Q . Again, we found the same problem encountered for HD 
14 947. The best fit of H 7 gives a mass-loss rate different from 
the best fit of H Q , without any possibility to fit both simulta- 
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Fig. 18. As Fig. 14, for HD 15 629. 
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neously, with a discrepancy of roughly a factor of two. H Q in- 
dicates a value of about log M = -4.75 and H 7 one around log 
M = -5.0, but this low value of log M would also imply that 
the gravity is lower than in the case of higher mass-loss rate. 
This is shown in Fig. 21, where we see the fit for H 7 at log M 
= -4.75 and -5.0 at log g = 3.80. Comparing with Fig. 19 we 
see that the effect of lowering the mass-loss rate is larger than 
any change of gravity at log M = -4.75. At the lower log M of 
-5.0 we would obtain a lower and more constrained gravity, 
but it is impossible to obtain even a moderate fit to H Q . In any 
case, the mass-loss rate is very large, although we should point 
out again that our mass-l oss rate coincides well with other H Q 
determinations (Puis et al. 1996; Lamers & Leitherer 1993) but 



it is la rger t han that derived from radio fluxes (Lamers & Lei- 
therer 1993 ) by a factor of 4 if we adopt the value from the fit 
of H Q , and by a factor of 2 if we take the one more consistent 
with H 7 . 

The temperature is now much lower than in the plane- 
parallel case, as for HD 14947, again due to the combined ef- 
fects of sphericity and the new criterion for the He ionization 
equilibrium. In agreement with this, the fit to Hell AA4541, 
4686 is poor. The formal change in log g is large with respect 
to the plane-parallel values (0.20 dex) and now HD 15 570 has 
also the largest spectroscopic mass, even formally exceeding 
100 M Q . Little can be said about the actual mass and the mass 
discrepancy. In Table 3 we see that the spectroscopic mass is 
30% larger than the evolutionary one, but we did not attemp to 
bring both into agreement (which could be possible adopting 
log g = 3.65), since for HD 15 570 gravity and mass are only 
formal values in the centre of a large uncertainty area, and have 
to be regarded as rather inaccurate. It is clear however that HD 
15 570 is one of the most extreme O stars in the Milky Way, 
and probably one of the most massive stars. 

5.7. CygOB2#7, 03 If 

The line fit for this star, displayed in Fig. 22, shows similar 
problems to those of HD 14 947 and HD 15 570 discussed be- 
fore. The mass-loss rate is quite large, again exceeding 10~ 5 
M©/yr, and the parameters point to a very massive star. The 
large mass-loss rate results in a change of log g with respect to 
the plane-parallel hydrostatic models, which is relatively large 
(0.15). However, the mass discrepancy is still present, since the 
spectroscopic mass is still a factor of 2 lower than the derived 
evolutionary mass. Another important change with respect to 
plane-parallel models is the new helium abunda nce, n ow e= 
0.18. It is interesting to note that Herrero et al. (1999) could 



not find evidence of any helium discrepancy in their analyses 
of other, less extreme Cyg OB2 stars. 



Fig. 19. The line fit to H 7 in HD 15 570 with log g = 4.05 (solid 
line) and log g = 3.55 (dashed line). All other parameters are as 
in Table 3. The abscissa gives wavelength in A and the ordinate 
relative fluxes. 



6. Additional considerations 

6.1. The H-y problem 

As we have seen in the foregoing analysis, for three of our 
objects with dense winds the wind emission seen in H 7 is in- 
consistent with the mass-loss rate derived from H a . It may be 
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Fig. 20. The line fit for HD 15 570. The fit shown corresponds 
to the model fitting the red wing of H Q (see text). 



questioned why this obvious problem has not already been dis- 
covered by Puis et al. (1996) who found no such discrepancy 
and derived — for the object common to both investigations 
(HD 14947Q) — a lower values for log g than follows from 
our analysis. 

This difference has the following origin. For the problema- 
tic objects, the Doppler core of H 7 becomes optically thick in 
the wind. Thus, it depends only on the wind conditions (mainly 
in the transonic region) and is independent of any underlying 
photospheric radiation. Consequently, it can be correctly fitted 
by the approach described in Puis et al. and would compare to 
our results. 

On the other hand, the wings of H 7 depend on the illumi- 
nating radiation from below the transonic region. In the concept 
applied by Puis et al., this radiation was taken from photosphe- 
ric profiles calculated on the basis of hydrostatic models, and 
the wind contamination was correctly accounted for. The fitted 
log g value was finally corrected for the difference between 
hydrodynamic and hydrostatic atmospheres, accounting for the 



: Note, that HD 15 570 has not been analysed by Puis et al. by 
means of detailed profile fits. Instead, only previously published va- 
lues for equivalent widths have been used, so that any comparison is 
obsolete. 



Fig. 21. The line fit to H 7 in HD 15 570 with log M = -4.75 
(solid line) and log M = -5.00 (dashed line) at log g = 3.80. All 
other parameters are as in Table 3. Comparing with Fig. [T^ we 
see that for this strong wind, changing the mass-loss rate has 
a larger effect than changing the gravity. The abscissa gives 
wavelength in A and the ordinate relative fluxes. 



different formation depths in an approximate and global way, 
where the correction turned out to be only moderate in the pa- 
rameter space considered. 

This procedure, however, is only justified if the major part 
of the wings is actually formed in some pseudo-hydrostatic en- 
vironment, which is the case if the wind densities are not too 
high. For increasing wind density, however, the difference be- 
tween hydrostatic and hydrodynamic stratification (e.g., Puis et 
al., Fig. 16) becomes increasingly larger, and the point where 
the transition between both regimes occurs shifts to correspon- 
dingly larger (mean) optical depths. Thus, for higher wind den- 
sities only the far wings are formed in a purely hydrostatic en- 
vironment, whereas the inner wings are severely affected by the 
wind conditions, which display a lower density at given opti- 
cal depth. Consequently, the use of hydrostatic profiles and the 
global correction applied by Puis et al. will inevitably fail under 
those conditions. Moreover, the derived values for log g will 
be too small compared to "reality", since, for a given log g, the 
hydrostatic densities are always larger than those in the wind 
regime. As a result, the deviations between the profiles calcu- 
lated by Puis et al. and ours should become largest just outside 
the Doppler core, and vanish at the extreme wings. Thus, it is 
not too surprising that the results of our consistent description 
deviate from the results given by Puis et al. with respect to gra- 
vity (derived from the wings), whereas the H 7 Doppler cores 
would be consistent with the mass-loss rate derived from H Q , 
however are unfortunately not visible due to rotational broad- 
ening. 
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Fig. 22. The line fit for Cyg OB2 #7 using a spherical model 
with the parameters given in Table 3. 



We like to point out that for a (small) number of objects 
with denser winds the determination of the stellar mass be- 
comes completely impossible, since in those cases even the 
continuum is formed in the wind and the reaction of any profile 
on pressure scale height becomes impossible. 

In order to check under which conditions this problem will 
arise, we have calculated the (minimum) continuum optical 
depth given by electron scattering as function of wind parame- 
ters: 



s e p(r)dr, 



0.4 cm 2 / g 



1 + IuoY 
1 + 4F 



(1) 



with s e the electron scattering mass absorption coefficient, Y 
the helium abundance (by number) with respect to hydrogen 
and Jhc the number of free electrons provided per He atom. By 
means of the equation of continuity and for the f3 velocity law 
which is used also in our model simulations, 



«(0 = «~(i-— f, b = i-(v / VQO y/p 



(2) 



(where v » 0. lt> SO und is the velocity at the transition point be- 
tween pseudo-hydrostatic and wind region, cf Santolaya-Rey 



et al., |1997| , their Sect. 2.2), we have 

Ms, 



T e (v) 



1 



47ri?*v 00 bp J v/v 



(3) 



where v' = v(i)/voo. Inserting typical parameters and deno- 
ting by A the optical depth-like quantity 



.4 



M 



IORq 1000 km/s 1 + In e Y 



10- 6 M Q /yr R* Voo 1 + iY ' 
the electron scattering optical depth is finally given by 



T e {v) w 0.028 A (-In ) for/3 =1 

Voo 



(4) 



0.028- 



— - -(1- — — ), £V 



(5) 



(approximating b by unity). In Fig. ^3| we have plotted this 
quantity, evaluated at v/voo = 0.001 (which compares roughly 
with the velocity at the transition point for hot stars) as a func- 
tion of A and (3. The dotted line gives the optical depth of 2/3, 
and for all objects with r e (A) > 2/3 the continuum is defini- 
tely formed in the wind. We have indicated the position of our 
sample stars by asterisks, where the number refers to Table |[ 
Obviously, object no. 2 ( = HD 15 570) lies just at the border 
line, and thus a gravity determination is almost impossible (cf., 
Fig. 19 for the influence of gravity on H 7 ). For the objects nos. 
1 (Cyg OB2 #7), 5 (HD 14 947) and 6 (HD 210839) at least the 
inner wings of H 7 are severely affected by the wind, and thus 
result in a larger gravity than found by Puis et al. (if all other 
parameters remain the same). For the remaining three stars, the 
continuum is formed solely in the hydrostatic part, and the de- 
rived numbers should coincide with the approximate method, 
as is actually the case. 

6.2. Influence of Hell resonance lines 

Those readers in particular who are not well familiarized with 
spectral analyses of hot stars may question why we did not take 
the most prominent line of Hell in the optical part of the spec- 
trum, namely Hell 4686 (n = 3 — * 4), into consideration so far 
(although we have plotted it for all models and the discrepancy 
is obvious). This question is completely legitimate, since the 
upper level of this transition is just the lower one of our strate- 
gic lines He II AA 4200, 4541, and should be reproduced with 
a similar degree of precision if our models were reliable. 

However, it is well known that this line (if formed in the 
wind) is extremely difficult to fit, and, to our knowledge, has 
never been used in any kind of NLTE-analysis of luminous O- 
stars. Usually, if one compares the predicted profiles to obser- 
vations, the synthetic line turns out to be too strong in emission, 
even if all other lines including H Q do perfectly fit. 
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Fig. 23. Electron scattering optical depth at transition veloc- 
ity v/Voo = 0.001 as function of A (cf., Eq. g) and /3 = 1.3 
(dashed-dotted), 1.0 (fully drawn) and 0.7 (dashed). The posi- 
tions of our programme stars (according to the derived wind 
parameters) have been indicated by asterisks. 



This rather unsatisfactory behaviour, which is normally by- 
passed by simply excluding Hell 4686 from the line list, re- 
lates to the extreme sensitivity of the participating levels on 
the treatment of the He II resonance lines and their sensiti- 
vity to line-blocking (for a discussion concerning this problem 



of the formation of He 1-lines, cf., Santolaya-Rey et al. 1997; 
see also the related discussion concerning the ionization struc- 
ture of WRs by Schmutz 1997 1. In standard simulations for 
wind conditions as described here, where line-blocking effects 
are excluded, the dominating background "opacity" below 303 
A (referring to Hell Ly a ) is Thomson scattering, leading to 
extremely enhanced radiation temperatures at the resonance- 
line' frequencies. Compared, for example, to a detailed balance 
situation (see below), the ground-state becomes depopulated, 
which in consequence (and in connection with the increased 
escape-probabilities due to the velocity field) prohibits He III 
from recombining and gives rise to much weaker absorption 
edges at 229 A, c ompared to plane-parallel simulations (cf., 
Gabler et al 



1989 



, especially Fig. A2). 
Moreover, since the radiation temperature is increasing to- 
wards higher frequencies (due to the decreasing bf-opacity at 
lowest photospheric levels), the NLTE departure coefficients 
are larger for higher levels than for lower ones. Thus, in addi- 
tion to the wind emission by geometrical effects, the lines be- 
tween excited levels (predominately He II 4686) are contami- 
nated by a strong source function oc 6 U ppor/^iowcr > 1, which 
leads to a much stronger total emission than would be the case 
if the resonance lines were of less importance. Actually, a pi- 



Table 5. Different approximations for treatment of Hell res- 
onance lines: The value of k corresponds to the definition in 
Eq. (||), and the line styles are the same as in Figs. ^4| and 26. 
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"standard" 






dotted 
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detailed balance 






- dashed 
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blocked flux 
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1Q -14 


fv fully drawn 
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blocked flux 
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1Q -16 


fl dashed-dotted 



shown that the emission of He II 4686 could be significantly 
reduced if line-blocking was accounted for correctly. 

In order to investigate in how far the above effects are of 
influence for our analysis (especially for the strengths of He II 
4200 and 4541, respectively), we have run a number of simula- 
tions with different treatments of the He II resonance lines, for 
the example of our final model for HD 14947, where the majo- 
rity of lines is formed in the wind. To check our hypothesis that 
the dominating effect leading to erroneous results follows from 
the increased pumping by resonance lines, two principally dif- 
ferent approaches were considered, which should give similar 
results if the hypothesis were correct. 

On the one hand, we set all Hell resonance lines into de- 
tailed balance. Alternatively, we simulated an additional p- 
square dependent background opacity in the decisive frequency 
range 227 A < A < 400 A, defined by 



add 



cP ■ fv, 



(6) 



with different values of n between 3 • 10 -14 and 1 • 10~ 16 . The 
frequential dependence /„ was assumed to be either increasing 
or decreasing, 



f 1 

J V 

q» = 



i - q v 
0.5(14 
l/A- 



1/227 



1/400 - 1/227' 



(7) 



lot investigation by Sellmaier (1996) for the case of £ Pup has 



and the appropriate emission component has been set to Planck. 
The value of n corresponds to the inverse of the electron den- 
sity at that point where the additional opacity reaches the same 
value as the electron scattering opacity. From the numbers 
given above, it is obvious that our choice is rather low com- 
pared to what might be expected in reality. For our final dis- 
cussion, we have selected four models with parameters given 
in Table pi 

Fig. e3 verifies the expected behaviour for the Hell ground- 
state. Both for the models with detailed balance as well as with 
simulated background opacities, Hem begins to recombine in 
the outer atmosphere as long as the strong upward rates present 
in model 1 (dotted) are no longer active, so that a significantly 
enhanced ionization edge develops. Only for the model with 
the lowest value of k = 1 • 10~ 16 (not displayed), the influence 
of the background opacity becomes so weak that the model 
remains ionized throughout the wind. 
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Fig. 24. Radiation temperatures of emergent fluxes for the 
model of HD 14947, with different treatment of Hell resonance 
lines. Line styles as defined in Table B. 
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Besides the reaction of the ground-state, also the expected 
behaviour of the excited states (reversal of population) takes 
place as long as the radiation temperature does not significantly 
increase towards higher frequencies (model 3) or the resonance 
lines are not active (model 2). A comparison between the co- 
rresponding Hell 4686 profiles and the observations (Fig. |25|) 
shows that our simulations are in almost perfect agreement, 
whereas our standard model displays much too much emission. 
It is interesting to note here a remark given by our anonymous 
referee. Henrichs (1991) reports that the equivalent width of 
He II A4686 in ACep varies in concert with the high-velocity 
edge of the C IV A 1550 line, suggesting that both variations 
share a common origin that could rely on the behaviour of the 
background opacity or the resonance lines of He II, as studied 
here. 

A final check on how far our models compare to reality al- 
lows the comparison with the UV-line Hell 1640 (n = 2 — > 3), 
since the lower state of this transition behaves differently from 
the other excited levels. In those cases where the resonance 
lines are no longer active (detailed balance) or are of negli- 
gible importance (background opacities) and the Hell Lyman 
edge becomes optically thick throughout the atmosphere, this 
state becomes the effective ground state of the ion^[ Thus, its 
population is predominantly controlled by the photoionization 
balance at 911 A coinciding with the hydrogen Lyman edge. 
Since this edge is optically thin, the 2nd level becomes strongly 
overpopulated because of the diluted radiation field, i.e., the ab- 
sorption should be larger than for the standard model. A com- 
parison with the observed IUE profile^] shows that our models 
are on the right track, in contrast to the standard model which 
predicts too little absorption. Only at higher velocities, i.e., in 



4 This situation corresponds to the behaviour of H a in A-type su- 
«rgiants. 

5 SWP 10724 
and R.K. Prinja. 



pergiants. 

5 SWP 10724, kindly provided in reduced form by I.D. Howarth 



Fig. 25. Comparison of the He II 4686 profile observed in 
HD 14947 with those obtained in our simulation of line- 
blocking (full drawn), of detailed balance for He II resonance 
transitions (dashed) and in our standard model (dotted). 
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Fig. 26. IUE spectrum of HD 14947 5 , HellA1640, compared to 
profiles generated from models with different treatment of Hell 
resonance lines. Observations shifted by 90 km/s to the blue, 
velocity dispersion Uturb =100 km/s applied to account for the 
red-shifted emission peak (see text). Line styles as defined in 
Tabled. 



the outermost wind, are they too strong compared to observa- 
tions, whereas for the inner wind, which is the decisive part 
concerning our analysis, they are in perfect agreement. A com- 
parison with model 4 shows that some fine-tuning might im- 
prove even the situation at larger velocities. (We note that in 
order to fit the position of the emission peak, we had to apply 
an artificial velocity dispersion of roughly 100 km/s, consistent 
with the values found from the analysis of UV resonance lines.) 
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In conclusion, we found that suppressing the large upward 
rates from resonance lines present in our standard model gives 
rise to a different population of ground and excited levels. The 
exact mechanism for this suppression, however, seems to be 
irrelevant to the results, and the synthesized profiles of pre- 
viously problematic lines compare well with the observations. 

We are now able to check the consequences of the mani- 
pulations outlined with respect to the strategic lines analysed 
so far and to derive constraints on which lines are more robust 
concerning our present ignorance of the real situation. 

Fig. H shows the profiles of He II AA 4541, 4200 produced 
by the different simulations and the "standard" model. We see 
that the wings of the He II lines become stronger in the sim- 
ulations as a consequence of the now reduced departure coef- 
ficients of the upper levels. However, the effect is smaller in 
the He II A4200 line (transition 4 — > 1 1) than in the He II A4541 
(transition 4^9) because transitions involving higher levels 
are weaker and form closer to the photosphere, so that the 
increased source function is not so visible. Since He II 4200 
shows up to be more stable, this is the preferred line in case 
of any discrepancy. This is a change of criterion with respect to 
former analyses in our group, but we prefer always to follow a 
single criterion that allows us to understand physically changes 
in the derived parameters. In addition, we should mention that 
also H Q becomes weaker, especially in the blue wing, as a con- 
sequence of these effects in the overlapping He II line. That is 
the reason why we prefer to fit the red wing in case of diffi- 
culties like those in A Cep or HD 15 570 (and in concert with 
the findings by Puis et al. (1996 1, who had also to manipulate 
the He II departures predicted by unified models if the wind 
was strong). Finally, we point out that He I is not so strongly 
affected if these lines are formed purely in the photosphere, 
but, as has been shown by Santolaya-Rey et al. (1997) they 
are also influenced by any effects that modify the population of 
the He II ground level (see Sect. 3.4.1 in Santolaya-Rey et al., 
where this problem has already been discussed). 

7. Discussion 

The first thing which is evident from the spectral analysis is 
the increasing difficulty in fitting the spectra of the earliest 
types with plane-parallel model atmospheres. In particular, for 
50 000 K and above we were unable to determine the stellar pa- 
rameters in the same way as we did for stars in Paper I and for 
the rest of the objects in the present study. This is true even for 
the relatively large gravity star Cyg OB2 jfl, and was nearly 
the case for the relatively cool star HD 15 558. Plane-parallel 
analyses are still useful, however, because they apply up to 
temperatures of 50 000 K and because they can be used as con- 
straints for the analysis in the larger parameter space demanded 
by more sophisticated models. 

The temperature scale defined by plane-parallel, hydrosta- 
tic, non-blanketed models is probably too hot (see Vacca et al. 
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1996; Harries & Hilditch 1998). For example, Hubeny et al 



(|1998j) have shown that the same quality fit can be achieved for 
10 Lac with a line -blanketed model at 35 000 K and with the 



Fig. 27. Line profiles of He II 4541 (left) and He II 4200 (right) 
for the simulation of line-blocking (fully drawn), of detailed- 
balance (dashed), and for the "standard" model (dotted). Cal- 
culations were made for the final model of HD 14 947. 



non-blanketed model used in Paper I at 37 500 K (all other pa- 
rameters remaining the same). The effect can be even larger for 
the stars analysed here, since the line-blocking effects in plane- 
parallel, hydrostatic models move the stars towards higher tem- 
peratures. 

In the spherical, non-hydrostatic models we have used here, 
line-blocking has been simulated, and we have found that it has 
the same effect as to keep the He II resonance lines in detailed 
balance (which were already in detailed balance in the plane- 
parallel models, so that they did not show the influence on the 
He II lines of the Pickering series we observe in these spherical 
models). This effect had important consequences in our ana- 
lyses. It led us to change the temperature criterion and adopt 
He II A4200 as the main line for the fit, which resulted in lower 
temperatures than if we had adopted the Hell A4541 in those 
cases in which we could not fit both lines simultaneously. The 
simulations indicate that the reduction of the strong upward 
rates of the He II resonance lines (whatever the real physical 
cause) will help to bring both lines into agreement (although 
at the moment we cannot say whether it will bring them com- 
pletely into agreement). It also suggested to us that we should 
fit the red wing of H Q whenever the fit to the whole line was 
impossible. What could be important for the future is that we 
have shown that He II A4686 is also affected and its fit is highly 
improved (although still qualitatively) when the resonance lines 
of He II are kept in detailed balance, or the departure coeffi- 
cient of the second level is kept below its detailed balance value 
through additional background opacity. 

The plane-parallel spectroscopic masses are as usual lowe r 
than evolutionary masses (see Paper I or Vacca et al 



1996) 



The new temperature criterion does not strongly influence the 



A. Herrero, J. Puis & M.R. Villamariz: The Upper HR Diagram 



21 



mass discrepancy. In Paper I we showed that lowering the tem- 
perature (for whatever reason) will not bring the masses into 
agreement. In the present case, however, we find the mass dis- 
crepancy even for luminosity class V stars. This finding is con- 
trary to the conclusion we obtained in Paper I (and is in agree- 



ment with other authors; see, for example, Vacca et al. 1996) 



namely that luminosity class V stars do not show a mass dis- 
crepancy. The reason for this apparent contradiction is that the 
gravities we derive here are also low for these luminosity class 
V stars. Thus the conclusion should be rather that 'high-gravity 
stars do not show a mass discrepancy,' where the term 'high' 
actually depends on the strength of the radiation field. 

This indicates a problem of the hydrostatic models due 
to the intense radiation pressure, and in fact larger gravities 
are derived when using the spherical models with mass-loss, 
largely reducing the mass discrepancy, which is now about 
50%. This result is even reinforced when we realize that most 
temperatures are now lower due to sphericity and a new tempe- 
rature indicator. An additional contribution came in some cases 
from the fact that the wings of H 7 can be strongly affected by 
wind contamination. In one case (HD 15 570) this contamina- 
tion is so strong that actual information about the gravity from 
the wings of H 7 is lost. In all other cases, however, the syste- 
matic effect that spectroscopic masses (without line-blanketing 
or blocking) are lower than evolutionary ones (without mixing 
mechanisms) is still present (note that formally the error bars 
overlap, and it is only because the effect is systematic that we 
can speak of a mass discrepancy). 

The problem can be alternatively formulated using the es- 
cape velocities. (Table 6 gives the escape velocities and related 
values). Escape velocities have been derived using log g va- 
lues uncorrected for centrifugal forces, i.e., the measured va- 
lues implicitely including the centrifugal force acceleration (cf. 
Sect. Q) which determine the effective escape velocities. 

Fig. shows the correlation between escape velocities and 
terminal wind velocities. We see that the diagram using evolu- 
tionary masses shows a good linear correlation between both 
quantities, whereas that with the spectroscopic masses shows 
a weaker correlation. However, the last diagram is in better 
agreement with the theory of radiatively driven winds. This the- 
ory predicts for the O-star domain (i.e., if the force-multiplier 
parameter 5 is small, cf. Friend & Abbott 1986; Kudritzki et al. 
198^) 



2.24- 



1 



with v csc the escape velocity and a one of the line force mul- 
tiplier parameters (the coefficient in the exponent of the line 
strength distribution function). With typical values for OB stars 
of a= 0.6-0.7 we obtain values of 3.3-5.2 for the ratio of 
terminal-to-escape velocity. This range of values is in agree- 
ment with those found here for the spectroscopic masses (see 
Table ^ note that the only point deviating from this range co- 
rresponds to HD 15 570, whose gravity and spectroscopic mass 
are very uncertain, and that HD 5 689, which is the leftmost 
point in the upper part of Fig. ^8[ would lie in the middle of the 



range if the absolute magnitude of a luminosity class III object 
were assumed). 

On the other hand, we derive very low ratios between ter- 
minal and escape velocities when using evolutionary masses 
(see Table Q), with an average ratio of 2.57, corresponding to 
a= 0.53 (which would be very low for the considered spec- 
tral range). Thus, Fig. ^8] seems to indicate that the evolution- 
ary masses are systematically too large (via the corresponding 
escape velocities), whereas the spectroscopic ratios are closer 
to the theoretically expected range. In Fig. 29 we can see this 
res ult, a lready contained in f ormer papers (Groenewegen et 
al. 1989; Lamers & Leitherer 1993), from a slightly different 



point of view. Here we have plotted the ratio of evolutionary- 
to-spectroscopic mass versus the ratio of spectroscopic-to- 
evolutionary a. We see that, except for the odd case HD 5 689, 
there is a strong correlation, indicating the relation between 
mass discrepancy and our present knowledge of radiatively 
driven winds. (Note moreover that HD 5 689 would perfectly 
fit into the correlation when an absolute magnitude of -5.78 is 
assumed, as appropriate for an 06 III star, but note also that 
the terminal velocity for this star was derived from the spec- 
tral type-terminal velocity relation from Haser 1995, as was 



that of Cyg OB2=#=7). We should stress that we adopt a /3-law 
for the wind velocity and determine then the a's directly from 
the derived relation between P^and escape velocity. Thus the 
agreement of the spectroscopic a's with the predictions of the 
radiatively driven wind velocity results in mutual support. 

At present it is still unclear whether the inconsistency found 
is due to some physical effects that should be incorporated into 
the evolutionary or into the atmospheric models, although the 
former are perfectly able to explain the discrepancy, at least 
qualitatively, when introducing mixing effects (Heger 1998| ) 
that could also affect surface abundances. 

Spherical models with mass-loss do not contribute to re- 
duce the helium abundances found with p lane-parallel m odels . 
In previous analyses (Herrero et al. 1995 ; Israelian et al. 1999| ) 
in which we used He I AA4387, 4922 and He II A4541 as tem- 
perature indicators, the derived helium abundances were simi- 
lar. As a result of the new indicator, the helium abundances had 
now to increase. However, we see that the stars of larger gravity 
(HD 15 558 and HD 15 629) do not seem to show He overabun- 
dances. Normal supergiants lie between the former and the fast 
rotators, that show the largest overabundances, in agreement 
with results from Paper I and the possibility that rotation plays 
a fundamental role in the chemical evolution of single stars. 
However, it is not clear whether a difference exists in the mass 
discrepancy between the rapid rotators and the other stars, or 
a correlation between the mass and helium discrepancies, that 
could indicate an overluminous evolution, as predicted by evo- 



lutionary rotating models (see Langer & Heger 1998, Maeder 



1998 or Meynet 1998 and references therein). We should stress 



here that we do not consider HD 5 689 as really showing such 
a large mass discrepancy as it appears to do in Tables 2 and 3, 
but that it is a problem of the stellar classification (or a problem 
of assigning the star to the Cas OB7 association). 



22 



A. Herrero, J. Puis & M.R. Villamariz: The Upper HR Diagram 



Table 6. T values, terminal velocities and escape velocities (in km s x ) obtained using the spectroscopic (sp) and evolutionary 
(ev) masses. In both cases we have taken into account the effect of the centrifugal force in reducing the escape velocity. 



Star 




V(sp) 


T(ev) 




Vool V CBC (sp) 


v CBC (ev) 


Vool v CBC (ev) 


Cyg OB2 #7 


2900 


0.71 


0.36 


586 


4.95 


1231 


2.36 


HD 15 570 


2600 


0.31 


0.43 


1149 


2.26 


882 


2.95 


HD 15 629 


3000 


0.45 


0.28 


781 


3.84 


1148 


2.61 


HD 15 558 


2800 


0.46 


0.43 


981 


2.85 


1039 


2.69 


HD 14947 


2400 


0.32 


0.25 


789 


3.04 


943 


2.55 


HD 210 839 


2250 


0.33 


0.31 


734 


3.07 


778 


2.89 


HD 5 689 


2500 


0.40 


0.10 


405 


6.17 


1013 


2.47 

















5 










r 6000 
























I 

l 4000 














4 








% E 


A 2000 


- ' — I — 






* 
















8 n 

> 

( 














s 2 








% -. 


) 


500 


100C 
[km s _1 ] 




1500 


1 




+ 

m 

































15C 



[k 



m s 



Fig. 28. The escape velocities obtained from spectroscopic 
(above) and evolutionary (below) masses, against the wind ter- 
minal velocities, all in km s _1 . The lines have the slopes pre- 
dicted by theory for a values of 0.6 and 0.7, e.g., 3.3 and 5.2. 
Plus signs mark the position of HD 5 689 if this star is assigned 
a magnitude of -5.78, typical of a 06 III star. Typical error bars 
have been plotted in the upper right corner. The abscissa error 
of Cyg OB2#7 in the upper plot is twice the corresponding 
error bar, due to its high T(sp) value 



Derivation of the mass-loss rates has been hindered by our 
finding of the inconsistency between H a and H 7 in those cases 
in which the wind is particularly strong. Thus we have found 
that both lines can demand mass-loss rates that differ up to a 
factor of two. We prefer the values given by H Q because it is 
more sensitive to mass-loss, and attribute the problem to diffi- 
culties in describing the wind in the transition zone. IR obser- 



0.80 0.90 1.00 1.10 1.20 1.30 1.40 
a(sp)/ot(ev) 

Fig. 29. The ratio a(sp)/o;(ev) against the mass ratio M ov IM s . 
Asterisks refer to the values quoted in Table 3 (in particular, 
the asterisk in the upper right corner corresponds to HD 5 689), 
whereas the cross marks the position of HD 5 689 if this star is 
assigned a magnitude of -5.78, typical of an 06 III star. 



vations and analyses giving information about this zone should 
help in the future to solve the problem. 

From the three stars for which we had some information 
about mass-loss rates derived from radio emission, only one 
(for which in addition only an upper limit from radio fluxes is 
available) shows agreement between the H Q and radio values. 
For the other two, radio mass-loss rates are a factor of three to 
four lower than H a mass-loss rates (and thus H 7 values would 
be in between). However, several facts should be taken into ac- 
count before claiming that there is a contradiction between both 
sets of values. First, the concerned stars are only two particu- 
lar, rather special cases (one, HD 15 570, is a very extreme O 
star, and the other, HD 210 839, is a strong non-radial pulsator, 
rapid rotator with a probably non-spherically symmetric wind, 
cf. Sect. 5.2); secondly, what we actually obtain from the obser- 
vations are Q values, proportional to (M/(Rvoo) 3 / 2 ) for H Q , 
and to (M/(R 3 / 2 Voo)) for radio and thus the derived values 
depend on the particular set of chosen stellar parameters; and 
thirdly, also the proportionality constants in the case above de- 



pend on model details (for example, Lamers & Leitherer 1993 
assume that He is doubly ionized in the wind of HD 210 839 
between 10 and 100 stellar radii, whereas our calculations in- 
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Fig. 30. The wind momentum-luminosity relation for the stars 
in our sample. The abcissa is log(i/L Q ) and the ordinate is 
the logarithm of the modified wind momentum (MWM), log 
(M VqqR - 5 ). Asterisks represent OB supergiants from Puis et 



al. ( 1996 ). Open diamonds, squares and triangles are respec- 
tively OB supergiants, giants and dwarfs from the present work. 
The positions of HD 14 947 and HD 210 839 in both studies are 
joined by a line, where the major difference results from the dif- 
ferent effective temperatures and helium-abundances attributed 
to the objects. 

dicate that from 28 stellar radii upwards, He is single ionized, 
which would have an immediate effect on the derived radio 
mass-loss rates, increasing them). Thus, we cannot state that 
there is a general problem (or even a particular one), or decide 
which value we should prefer for each object. 

Having derived all parameters we can obtain the modified 
wind momenta of all stars for which the radiation-driven wind 
predicts a tight correlation with luminosity. Our results are 
given in Table 3 and plotted in Fig. together with the data 
given by Puis et al. (1996) for Galactic supergiants. We see that 
our points for supergiants agree well with theirs (we have dis- 
cussed the differences in the individual studies). We have also 
plotted the regression lines derived from the supergiants from 
Puis et al. and from all the supergiants (the former and ours) to 
stress this point. We see that both regression lines are quite pa- 
rallel, indicating that the new values do not significantly change 
the known WLR. 



8. Conclusions and future work 

We have presented the spectra of seven Galactic luminous 
O stars, that have been analysed by means of non-LTE, 
plane-parallel, hydrostatic model atmospheres, including line- 
blocking, and by means of spherical, mass-losing models. 
We have used these analyses to study some additional effects 



(physical as well as due to applied approximations) that have 
an influence on the results. 

We have shown that line-blocking has to be taken into 
account when analysing very hot stars, in plane-parallel as 
well as in spherical models with mass-loss. Line-blocking in- 
directly affects the spectral-line diagnostics through changes in 
the level populations of He. Although for the spherical models 
with mass-loss we only made simulations, this is clearly one 
direction for future model improvements. In particular, we ob- 
tained very interesting results concerning the He II A4686 line 
and the blue wing of H a (actually, concerning the correspon- 
ding He II blend). 

One of the conclusions of the present paper is that around 
50 000 K (the exact value depending on gravity) we reach the 
limit of applicability of plane-parallel, hydrostatic models to 
real massive stars. This conclusion can be further refined, if 
we take into account the fact that line-blanketed models could 
give lower temperatures, or if we consider that spherical mo- 
dels with mass-loss will demand larger gravities. The final re- 
sult would be models less affected by radiation pressure, larger 
masses and a cooler temperature scale. 

Spherical models with mass-loss are thus needed to analyse 
these stars. We had to correct the gravities of all stars (except 
HD 5 689) when using these models, sometimes with incre- 
ments as large as 0.25 dex. Part of the differences in the ste- 
llar parameters derived from both sets of models are due to a 
change in the preferred temperature indicator. 

Our spherical models with mass-loss are also still not free 
from internal inconsistencies. For some models with strong 
winds, we see that HellAA4200, 4541 give different tempe- 
ratures, although we expect that inclusion of line-blocking will 
reduce this discrepancy, and will strongly improve the fits of 
He II A4686 and the blue wing of H Q . We have also found that 
for stars with strong winds there is a difference in the mass-loss 
rates derived from H Q and H 7 , that can reach a factor of two. 
We do not find good agreement with the two cases for which 
we have mass-loss rates from radio fluxes, although it is not 
possible to derive any firm conclusion from this fact. 

The helium and mass discrepancies found here are in agree- 
ment with the findings of Paper I. The only difference is that we 
find a mass discrepancy also for luminosity class V stars, but 
this is actually misleading. What we find is that the luminosity 
class V stars analysed here have gravities lower than usual for 
their spectral classification. The conclusion of Paper I should 
then be changed to state that we do not find the mass discre- 
pancy for stars with high gravities. At least a part of this effect 
is attributed to the fact that the intense radiation field affects 
the wings of the Balmer lines also for these stars. An additional 
possibility is that they actually neve r reac h the zero-age main 
sequence (see, for example, Hanson 1998 ). 

Spherical models with mass-loss largely reduce the mass 
discrepancy. Without solving it completely, the problem lies 
now in the systematic trend of spectroscopic values to be lower 
than evolutionary ones, but most of the individual values now 
agree within the formal error bars. We have also shown that the 
spectroscopic masses agree better with the predictions from the 
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radiatively driven wind theory, because they give ratios of the 
terminal wind velocity to the escape velocity (or equivalently, 
values of the a line force parameter) that are in the range pre- 
dicted by that theory. 

From the point of view of individual stars, we have ana- 
lysed some of the most massive and luminous stars in the 
Milky Way. We have found that three of them (Cyg OB2 
#7, HD 15 570 and HD 15 558) have particularly large initial 
masses, around or in excess of 100 M Q (depending on the tech- 
nique used for the mass derivation). Cyg OB2=^7, HD 14947 
and HD 210 839 have so extreme mass-loss rates that the wings 
of H 7 are strongly affected, and in HD 15 570 the wind is so 
strong that the exercise of deriving the gravity from the wings 
of H 7 results in highly uncertain values. On the other hand, 
we find that HD 5 689, the less luminous object in our sample, 
could have been wrongly assigned to Cas OB7, and might be a 
runaway star. 
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